function [x,y]=gampdfcheck(para1,para2);

close all;

b = para2^2/para1;
a = para1/b;

gaminv(.9999,a,b);

x = linspace(0,ans,100);
y = gampdf(x,a,b);

fmt = '%0.5g';
plot(x,y);
title(['Gamma with (' num2str(a,fmt) ', ' num2str(b,fmt) ')']);
xlabel(['\mu = ' num2str(para1,fmt)  '    \sigma = ' num2str(para2,fmt) ]);
